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Abstract 

We compute quark-antiquark pair production in the context of the Color 
Glass Condensate model for central heavy-ion collisions. The calculation 
is performed analytically to leading order in the density of hard sources 
present in the projectiles, and is applicable to quarks with a mass large 
compared to the saturation momentum. The formulas derived in this pa- 
per are compared to expressions derived in the framework of collinearly 
factorized perturbative QCD and in k±_ factorization models. We com- 
ment on the breaking of k± factorization which occurs beyond leading 
order in our approach. 



1 Introduction 

Heavy Quark pair production in the collinear factorization formalism of per- 
turbative QCD is a well developed subject [1, 2, 3], It is however not evident 
that this formalism is applicable in the kinematic regime where the heavy quark 
mass m is much smaller than the center of mass energy (m <C y/s). In this 
kinematic regime, small x effects may be important. An alternative formalism, 
the k± factorization formalism, was developed to describe the physics in this 
region [4, 5]. A key feature of this formalism is that, unlike the collinear factor- 
ization formalism, the small x gluons that produce the heavy quark pair have 
intrinsic transverse momenta on the order of the hard scale of the process. At 
high fcxj the collinear factorization formalism is recovered. A comparison of the 
fcj_ factorized formalism to recent hadroproduction data has been performed in 
Ref. [6]. 
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In Ref. [7], it was suggested that at small x, the intrinsic transverse mo- 
menta in the k± factorized formalism are of the order of the saturation scale 
Q s {x) [8, 9]. Since the saturation scale grows with energy, so does the intrinsic 
transverse momentum of gluons, thereby significantly increasing the hadropro- 
duction cross-section. The saturation inspired k± factorized formalism and phe- 
nomenological applications including comparisons with the NLO perturbative 
QCD formalism are reviewed in Ref. [10]. A related approach is the color dipole 
approach which has also been applied to heavy quark hadroproduction [11, 12]. 

A systematic way to study high parton density effects in QCD is the Color 
Glass Condensate (CGC) [13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 
27, 28, 29, 30]. Heavy quark production from the CGC was first discussed in 
the context of deep inelastic scattering [31] and subsequently in the context 
of inclusive and diffractive photoproduction [32, 33]. We will discuss here the 
problem of hadro-production of heavy quark pairs at very high energies. We 
will discuss specifically heavy quark production in the scattering of two large 
nuclei which has some simplifying features. Inclusive gluon production in the 
scattering of two nuclei has been studied previously, first to lowest order in a s 
and lowest order in the parton density [34, 35, 36, 37] and later to all orders in 
the parton density (and lowest order in a s ) [38, 39, 40, 41]. 

The results for inclusive gluon production in the CGC approach, at lowest 
order in the parton density, can be expressed in a kj_ factorized form [36, 37]. 
This k± factorization is however broken when inclusive gluon production is com- 
puted to all orders in the parton density [38, 39, 40, 41]. We expect a similar 
breakdown of k± factorization in heavy quark production at the small trans- 
verse momenta Aqcd < k± < Q s when parton density effects are computed to 
all orders [42]. For a recent discussion of k± factorization and the CGC, see 
Ref. [43]. 

We do not explicitly include the effects of quantum evolution in our ap- 
proach - the energy dependence of our results is determined entirely by the x 
dependence of the saturation scale Q s . An interesting consequence of quan- 
tum evolution in the CGC is the geometrical scaling [44] of distributions with 
Q s . This persists in the kinematic window Q s < k± < Q^/Aqcd outside the 
saturation region k± < Q s [45, 46, 47, 48]. A number of phenomcnological 
models have argued that data at HERA for structure functions, vector mesons 
and even charm [44, 49, 50, 51] (as well as data from RHIC [52] and NMC [53]) 
exhibit geometrical scaling. The competition between k± factorization preserv- 
ing quantum evolution effects and multiple-scattering kj_ factorization breaking 
effects immediately outside the saturation regime is a very interesting prob- 
lem [54, 55, 56, 57, 58] but we will not address it further here. 

This paper is organized as follows. We begin in section 2.1 by discussing 
the general formalism for pair production in a time-dependent external field. 
The pair production probability for a single pair is expressed in terms of the 
classical background field produced by the color charge sources in the two nuclei. 
The cross-section and the average number of pairs are obtained by averaging 
the pair production probability with the appropriate weight functional which 
governs the likelihood of different orientations of the color charge density. In 
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section 2.2, we argue that the leading order contribution to pair production 
is of order 0{p\p\), where p\ and pi are the color charge densities of the two 
sources. We write down explicit expressions, in covariant gauge, for the classical 
field that contributes to the probability at this order. These covariant gauge 
expressions were first derived by Kovchegov and Rischke [36]. For completeness, 
we reproduce their derivation in appendix A. With these expressions for the 
classical field, we write down an expression for the pair production amplitude in 
section 2.3. In the following sub-section, we show that this amplitude satisfies a 
Ward identity. The pair production probability is computed in section 3. Some 
details of the computation are given in appendix B. In section 4, we relate the 
results of our computation to the k± factorization approach. The breakdown 
of k± factorization is discussed in the following section. We conclude with a 
summary of our results. 

2 Pair production amplitude 

2.1 Basics 

In order to compute the pair production amplitude in the Color Glass Con- 
densate model, we first recall some results [59] about particle production in an 
external time-dependent 1 classical field. We denote G F (x,y) the Feynman (aka 
time-ordered) propagator of a quark from the point y to the point 2 x in the 
presence of the external color field, and G a (x,y) the retarded propagator of 
the quark between the same points. Both are Green's functions of the Dirac 
operator: 

(0 X - <d(x) - m)G F:R (x, y) = (x - y) , (1) 

but with different boundary conditions. For instance, the boundary condition 
for the retarded propagator is simply: 

Urn G R (x,y) = S(x-y) 7 . (2) 

The boundary conditions for the Feynman propagator are awkward and not 
very illuminating. 

From these propagators, it is customary to first perform a Fourier transform: 

G FiR (q,p) = J d 4 xd 4 y e^e-^G^y) , (3) 

and then to extract the "scattering matrix" Tf.r vi a: 

G F J q ,p) = (2nr5( q -p) G FR (p) + G° FR (q)T FR (q,p)G FR (p) , (4) 

1 More exactly, one needs that the external field be time dependent in any Lorcntz frame. 
For instance, a single moving nucleus, which generates a time dependent classical color field 
in the laboratory frame, cannot produce pairs. 

2 The order of the points may seem unnatural, but it is chosen so that the Dirac matrices 
(which should be read from the endpoint to the starting point of the propagator) appear with 
the correct ordering. 
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where G° is the free Feynman or retarded quark propagator. This definition 
removes the free term, and amputates the external legs of the propagator. 

Several quantities related to pair production in the external classical field can 
then be expressed simply in terms of these scattering matrices. The probability 
for producing a single qq pair is given by 3 : 

ft m = |(0 in |0 out )| 2 / / j^j- W q )T Aq ,- p)v ( P )f , (5) 

where the argument A 11 has been used to remind the reader that this is a quantity 
defined for one particular configuration of the external classical field. The field 
is a functional of the hard sources p\ and p 2 (one for each nucleus) that 
generate this classical field. The prefactor 4 |(0i n |0 O ut)| is the square of the 
overlap between the vacua at x = — oo and x = +00. This quantity is strictly 
smaller than unity in an external field that can produce pairs. It is necessary 
for unitarity to be preserved (for instance to ensure that the probability Pi [A^] 
is always smaller than unity). The cross-section for single pair production can 
be obtained from here by: 

<ri = J d 2 b J [Dp 1 Vp 2 ]W[p 1 ,p 2 ;b}P 1 [A»[ Pl ,p 2 ]] , (6) 

where W[pi, p 2 ]b] is the functional weight that defines the statistical distribu- 
tion of the hard sources at a given energy. This functional depends on the 
impact parameter b between the two nuclei. The latter must be integrated out 
in order to convert the probability P\ [A^] into a cross-section. 

Similarly, the average number of produced pairs in a given external field 
configuration is given by: 

The average number of pairs produced in a collision at a given impact parameter 
b is given by: 

n(b)= J[Vp 1 Vp 2 }W[p 1 ,p 2 ;b}n[A^[p 1 ,p 2 }} . (8) 

Note also that both Pi and n can be expressed as differential distributions by 
undoing the integrations over the momenta q and p of the quark and/or the 
antiquark. 

3 From eqs. (4) and (5), and given the fact that the canonical dimension of the spinors u(q) 
and v(p) is (momentum) 1 / 2 , one can check that Pi is indeed dimcnsionless. 

4 This prefactor is also the probability that no pair is produced during the collision, which 
explains why it has to be smaller than 1. Diagrammatically, this quantity is the sum of 
the "vacuum-vacuum" diagrams, that have no external legs besides those that connect to 
the external field. More precisely, if we denote V the sum of all connected vacuum-vacuum 
diagrams, then | (0i n |0 O ut) | 2 = exp(— 21m V) [59, 32, 33]. The expansion of ImV in powers of 
the hard color sources starts at the order 0(p\p^)- 
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When we give explicit expressions for the average over the distribution of the 
hard color source charge densities, p\ and p 2 , we use the McLerran-Venugopalan 
Gaussian model for the functional W[pi, p 2 ;b], which reads: 



W[p 1 ,p 2 ;b] =exp- / d 2 x± 



pl,a(x±)p lj g(x_ l _) p2,a( X ±)P2A X ±) 

2pl(x ± ) 2p 2 2 (x ± -b) 



(9) 



where the functions p\ (x±) and p 2 (x±) correspond to the color charges squared 
per unit area of the two nuclei. They are functions of the transverse coordinate 
that describe the number of hard color sources per unit area at the transverse 
coordinate x±. (The origin of the transverse coordinates is taken at the center 
of the first nucleus, hence the argument x± — b for the function that describes 
the second nucleus.) The canonical dimension of p\ 2 is (momentum) 2 . In a 
model of a nucleus with sharp edges, these functions would be almost flat inside 
the nuclear disc and zero outside. 

In the MV model, p 2 is energy (or x) independent. The x dependence comes 
in through quantum evolution [16, 17, 18, 19, 20, 21, 22, 24, 26, 27]. In the MV 
model, one can relate p? to the saturation scale Q s 5 : 



The relation between the two scales is more non-trivial under quantum evolution 
since the color screening scale is itself of order Q s [45, 46, 47, 48, 61]. We will see 
later in Section 3 that, in our formalism, the x dependence of pair-production 
enters only through the x dependence of the saturation scale. 

2.2 Leading order approximation 

The relations of eqs. (5) and (7) are completely general and are valid to all orders 
in the sources p\ and p 2 . In QCD, the external field A^[p\,p 2 ] is a nonlinear 
functional of the sources p\ and p 2 , that receives contributions to all orders in 
Pi and p 2 ■ In addition, the scattering matrices T F R are themselves functionals 
of the external field A 11 , receiving contributions to all orders in A^ . These facts 
make the quantities Pi [A 11 ] and n [A^] extremely complicated functionals of p\ 
and p 2 that cannot be calculated analytically. A complete solution would likely 
require a numerical solution [42]. 

The leading term in the hard sources can however be calculated in closed 
form. This approximation is justified if there is a scale in the problem (the 
quark mass or the quark transverse momentum) which is large compared to the 
scale set by the density of hard color sources (namely, the saturation momentum 
Q s ). We will assume that the strong coupling constant a s is sufficiently small 
to justify keeping only the lowest order in a s . However, all corrections of order 
a s ln(s) are included since they are resummed via the quantum evolution of the 
functional W[pi,p 2 ; b]. 

5 For a discussion of different conventions for /i 2 , see Ref. [60]. 



Q 2 s (x,x±) = a s N c p 2 (x,xj_)\n 




) 



(10) 
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One can convince oneself that both Pi [A^\ and n [A 11 ] start at the order 
0{p\p\) m tri e hard sources. This is so because these quantities are squares, 
and because the quark line must be attached at least once to each nucleus if the 
pair is to be produced on-shell. At this order in the hard color sources, there 
are several important simplifications that are worth mentioning. Firstly, the 
prefactor |(0i n |0 ou t)| 2 that appears in Pi can be replaced by 1. Indeed, one has 
I (0in|0 O ut) 1 2 = 1 + C(/°i/°2)- Keeping higher order corrections in this prefactor 
would only produce terms of order higher than 0{p\p\) in Pi. The second 
simplification occurs because there is no difference between the time-ordered 
and retarded amplitudes. Indeed, since the probabilities P2, P3, ■ ■ ■ to produce 
two or more pairs in a collision are of higher order, we find 

n[A"] = Pi [A»] + 2P 2 [A»] + 3P 3 [A"] + • • • 

w Pi [A»] . (11) 

The time-ordered and retarded amplitudes are the same at this order because, 
in the language of Feynman diagrams, none of the intermediate state quark 
propagators can be on their mass-shell (namely, the it prescription of the quark 
propagators does not matter). 




Figure 1: The leading contributions to the pair production amplitude in terms 
of A^, A% and A^ 2 . The gluon line terminated by a cross represents an insertion 
—igA^(x) of the external field. 

Our first task is therefore to obtain the classical field A^ [pi , p 2 ] up to the 
order 0{pip 2 ). We can formally write: 

A»(x) = A$(x) + A$(x) + A&(x) + ■ ■ ■ , (12) 

where A^(x) is the contribution of order 0{pi) to the classical field, A^ (x) the 
term of order 0(p2), and A^ 2 (x) the term of order 0(pip2)- In terms of these 
contributions to the classical field, we have to evaluate the three terms in figure 
1. The derivation of the classical color field in the covariant gauge up to the 
order 0(pip 2 ) can be found in [36]. We merely quote the results in momentum 
space here, and refer the reader to the original literature or to appendix A of 
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the present paper where this derivation has been reproduced in our notation 6 : 
A+Jk) = 2ng5(k-)^ Pl , a (k±) , A^ a (k) = A\Jk) = , 
A^ a (k) = 2ng6(k+)-^-p 2 , a (k ± ) , A+Jk) = A\Jk) = , (13) 

K I 



A + 12 Jk) = -igi^- J ^Ljfc+ + -|j A+ b ( kl )A 2 Jk 2 ) > 

A- 12 , a {k) = ig ^- J ^ |fc 2 - + if} A+ 6 (fci)^c(fe) , 

4 2 , a (*0 - iff^ J J^{k l 2-k\} A+JkJA^) , (14) 
where k 2 = k — k\. 

2.3 Pair production amplitude 

Following figure 1, we can write the leading order pair production amplitude as: 
u{q)T FR v{p) = M(q,p) = Mi +2 (q,p) + Mu(q,p) , (15) 
where we have defined 

M 1+2 (q,p)^ J ^(-igAlMX-igAZtfo)) 

xu(q) {lnt a G°{q - fci)7„t 6 + ^ v tbG°(ki - p)7„t„} v(p) (16) 

and 

M 12 {q,p) = (-igA^ a (p + q))u(q) lt ,t a v(p) , (17) 

where the t a ^ are color matrices in the fundamental representation. G°(k) is 
the free quark propagator 7 : 

G°(k) = i ,J + ™. , (18) 

and = p + q — k\ . 



6 The light-cone coordinates are denned as x^ = (x° ± x 3 )/\/2. With this convention, the 
invariant scalar product of two 4- vectors is a ■ b = a+6~ + a~fe+ — ax • bx and the element 
of 4- volume is d 4 x = dx + dx~ d 2 x±. 

7 The prefactor i in this propagator is purely conventional. However, it cannot be chosen 
independently of the —i in the field insertion —igA 11 . 
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Using the fact that the momenta p and q are on-shell, that all but one 
component of and are zero, and finally fcf = k^ = 0, we can obtain 
easily: 

74; 



M 1+2 (q,p) = tg 2 J ^^A+Jk^A^) 



, . f m — 7 i • (q I — fci_i_) i 

XU (q) 1~ - _ — f K \ r^f 2 7 + ^fc 
L 2<? p+ + (q ± - kij_) z + m z 

(19) 

We readily see that the denominators in this expression are positive definite, 
which justifies a posteriori the fact that the ie prescription in the free quark 
propagator is irrelevant at this order. Using eqs. (14), we can similarly write 
the term Mi2(q,p) as: 

-lit 



MMq,p) = i 9 >-^ J ^A+MA-^) 

xu(q)S^- (V - ^L) - 7 + (k- M^) + 7± ■ (fc 1± - fe 2± ) ji>(p) • 

(20) 

Here again, all the denominators are strictly positive, making the ie prescriptions 
irrelevant. Note that we can write the expression of A4u as follows: 

M 12 (q,p) = ig 2 -^^ J ^La+MA-^uMMp) , (21) 



{c + =p + + q +-J^-,c- = -^--p-- q -,c ± = k 2± -k 1± ) . (22) 



where d 1 is the following 4- vector: 

p~+q~ ' ' + 

We note that the vector c M is indeed nothing else but the standard effective 
Lipatov vertex (see eq. (80) of [62] for instance). In particular, one may check 
that it obeys the transversality relation: 

(p + q)-c = 0. (23) 

At this stage, it is trivial to replace A^ and A^ by their expression in terms 
of the source densities (eqs. (13)) in order to obtain the leading order pair 
production amplitude as a functional of the hard color sources p\ and p 2 ■ 

2.4 Ward identity verification 

One can write the pair production amplitude M(q,p) as follows: 



M(g,p) = J ^iyfe)^(feKt(fa,fe;9,p). 



(24) 
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The expressions of the previous subsection tell us that: 

m ah + ( fc i' fc 2;9,p) = ig 2 

xu(q) 



1 ( m — Y± ■ (q±_ - k!±))-f + t a t b 7+(m+7_L • (p ± - feu.))7 t b t a 



+ - 



2q p+ + (q ± - k 1± ) 2 + m 2 2q+p + (p ± - k 1± ) 2 + m 2 

[t a , tb] 



h 2 

^2 . 



h 2 

2 fc+ 



+ 7± ' (fel_L-fc2J_) 



U(P) , 

(25) 



(p + q) 2 

while the other components of m^X do not appear in the expression of the pair 
production amplitude. One notes that rn^(ki, k 2 ; q,p) is not the amplitude 
for the process gg — > gg. Indeed, contracting the 55 — > gg amplitude with the 
classical fields A^(ki) and A 2 (k 2 ) would only give the first three diagrams of 
figure 2, and leave out the bremsstrahlung diagrams. 

It is interesting to investigate the limit ku_ = of m~ b + (ki, k 2 ; q,p). We 
have: 

lim m~ b + (ki,k 2 ;q,p) = 

■ 2-/ J 7"("1-7_L -Q±h + tatb , 7 + (™ + 7_L ' P± )7"*6*a 

= «5 w(9)s — ^ — — — o 1- 



2q~p+ + q^+m 2 



2q+p~ + p 2 ± + m 2 



[ta, tb] 

\p + qf 



7 k~l — 7 4 



-7_L ' fc 2± 



>v(p). (26) 



Using the Dirac equations obeyed by the spinors u(q) and v(p), as well as the 
following relations: 

q\ + m 2 = 2q+q- , p\+m 2 = 2p+p- , 

kf =p + + q + , k 2 =p~ +q~ , k 2 ± =P± + q± (if k 1± = 0) 

7+7-7+ = 2 7 + , (27) 

we can simplify the above limit into: 



lim m ab + (ki,k 2 - 1 q,p) 



y+v(p) [ l a ' tb l ( 
' + q+ \ 



= ig 2 u(q)^ 
= . 

Similarly, one can check that: 



1 + 



p+ + q+ 



2(p-+q-) 



p + + q + 



lim m ab + (k 1 ,k 2 ;q,p) = . 



(28) 
(29) 
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These simple limits are consequences of the following Ward identities 8 : 

ki f ,m^(k 1 ,k 2 ;q,p) = k 2u m t 2(k u k 2 \q,p) = . (30) 
Indeed, since fcf = 0, the first one implies for instance: 

ktm- b + (k 1 ,k 2 ;q,p) = kimi+(k 1 ,k 2 ;q,p) = 0(k\) . (31) 

In other words, m~ b + {k\, k 2 ; q,p) must vanish linearly with ki± when ki± — > 
in order to fulfill the Ward identity. The physical meaning of these vanish- 
ing limits is that in order to produce the quark-antiquark pair on-shell, some 
transverse momentum has to come from both nuclei. As we shall see later, this 
property is also essential in order to soften some collinear singularities. 

3 Pair production probability 

3.1 Average over the hard color sources 

From eq. (24), we can write the pair production amplitude in terms of the hard 
color sources: 

M(q,p)=g 2 J (2?r)2 fc2 — 2 m ab + (k u k 2 ;q,p) . (32) 

Squaring this amplitude, and averaging over the sources, we obtain: 

™ = * m = 94 / wk; I A, I [VnVK] w WkM 

r d 2 k 1± d 2 k' 1± Pl,q(fcl±) p 2 .b{ k 2±) P*l.a>( k 'l±_) P*2.b'( k 2±) 
X J (27T)2 (2n) 2 k f ± fc2 ± fc* k f ± 

xTr (m-+(k 1 ,k 2 ;q,p)m-+*(k' 1 ,k 2 ;q,p)) , (33) 

where Tr denotes a trace of the color and Dirac matrices. In the McLerran- 
Vcnugopalan model where the functional W[pi, p 2 ;b] is given by eq. (9), the 
source averages are given in coordinate space by: 

J [V Pl Vp 2 ] W[p 1 ,p 2 ;b]p lta (x±)p 1 , a >(x' ± )p 2 , b (y ± )p 2t v(y' ± ) = 
= 6 aa '8 bb 't4( x ± )n 2 2 (y ± - b)8(x ± - x' ± )5{y ± - y'J . (34) 

8 Bccause of eq. (24), the object (fci, q, p) may be seen as the lowest order value of 
the correlator (Jf J^ipi/i) between the currents of the individual nuclei and a pair of fermions. 
The currents and behave like Abelian currents at the order we are considering (the non 
Abelian term in the current conservation law is of higher order in the sources pi^), which is 
the reason why we have Abelian-like Ward identities for m^ b {k\ , k^; q, p). 



10 



The Fourier transform of eq. (34) reads: 

j [Vp{Dp 2 ] W [p 1: p 2 ;b] Pl . a (k 1± )p 2 ,b( k 2±)P* 1 ,a>(k' 1 ±)p*2, b >(k2 ± ) = 

= 5 aa ' '5 W 'e- i6 -< fc ^-*i^(fci ± - feW)Mi(fei± - fcix) , (35) 
where /z^ 2 (fc±) is the Fourier transform of p,\ 2 {x±_). We have used the fact that 

= p± + <?± - fcix- 

3.2 Impact parameter and energy dependence 

At this stage, we have the following expression for the pair production proba- 
bility: 

P,(b)=n(b) = a* f ^ [ d3<1 f d2kl± ggLk e -*-(»ix-fcix) 
ll J ( J 9 J (27r)32c pJ / (2^2cJ (2tt)2 (2^ 

x^fci, fc' WW fr. , ( W a6 + ( fc i^2;g,p)m a -+ *(A^;g,p)) 

XMll K lX - K 1±JAM K 1-L ~ K 1^J ,2 7/2 ,2 ,/2 ' 

' C 1X' C 1X' C 2±' C 2X 

(36) 

For large nuclei, these Fourier transforms are strongly peaked around k±_ = 0, 
with a typical width of the order of 1/R where R is the radius of the nucleus. 
For instance, for R = 6 fm, the spread in transverse momentum would be of 
order Ak± <~ 1/R ~ 30 MeV. The factors n\ fc + contain a large momentum 
scale provided by the mass m of the heavy quarks, which we assume to be much 
larger than 1/R. Therefore, we can safely neglect the difference ki± — k' 1± 
in these factors. Regarding the denominators k 1A _k 1± k 2± k 2± , the situation is 
more delicate. Indeed, they could be arbitrarily small because of the presence 
of a collinear singularity when k\± — > or k\± — > p_|_ + q ± . We are going 
to assume that these collinear singularities are regularised by physics at the 
scale of the saturation momentum Q s 3> 1/R, so that we can also neglect the 
difference fci_i_ — fc'^ in the denominators. Therefore, we can easily integrate 
out the difference fcij_ — k 1± in order to get the impact parameter dependence 
of the pair production probability: 

Pi (6) = n(b) = g A j d 2 x^l( Xl _)pl{ Xl _ - b) 

f d 3 p f d 3 q f d 2 k 1± Tr (hab + (ki,k2;q,p)\ 2 ) 
X J {2irf2u p J {2nf2iu q j (2tt) 2 k\jz\^ ' ( ' 

The prefactor between the square brackets is the overlap of the density functions 
that describe the two colliding nuclei. In a crude model of "cylindrical" nuclei, 
it would be equal to p\p 2 S(b) where 5(6) is the area of the overlap between the 
two nuclei. The explicit expression of the trace that appears in this formula is 
given in the appendix B. 
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Since there is no quantum evolution in the MV model, the functions n\(x±) 
and /J,2(x±) do not depend on the longitudinal momentum fraction of the gluons. 
As a consequence, the differential pair production probability is boost invariant, 
in the sense that it depends only on the rapidity difference y p — y q of the quark 
and the antiquark. Strictly speaking, eq. (37) gives an infinite result because 
it contains also an integration over the mean rapidity (y p + y q )/2: this infinity 
should be cutoff by the fact that the rapidity of the produced quark and an- 
tiquark are limited by the rapidity of the projectiles. More formally, quantum 
evolution of the functional W[pi, p2',b] would introduce a dependence of the 
functions p1(x±) and p^(x±) on the gluon momentum fractions, such that they 
vanish when we reach the fragmentation region of the corresponding projectile. 

4 Relation to /c^-factorization and Collinear fac- 
torization 

In fcj_-factorized perturbation theory, it is well known that the leading order 
diagrams contributing to the production of a quark- antiquark pair are those 
listed in figure 2. The first two diagrams are the analogs of the diagrams one 
would have in QED, except for the fact that the coupling of the gluon to the 
quark line involves a non-commuting SU(3) matrix. The third diagram involves 
the 3-gluon vertex, and the fourth diagram 9 is a bremsstrahlung contribution 
required in order to ensure gauge invariance. 




Figure 2: The Feynman diagrams that can contribute to pair production at the 
order 0{p\p2). The bold lines represent the classical sources pi and pi. 

It is well known that the first two terms of figure 1 correspond to the first 
two Feynman diagrams of figure 2 in the limit where the transverse momentum 
of the two gluons is small in front of the longitudinal momentum of the hard 
sources. Moreover, it has been shown by Kovchegov and Rischke [36] that the 
term proportional to A% 2 (x) contains both the diagram with the 3-gluon vertex 
and the bremsstrahlung diagrams. Therefore, it is very interesting to compare 
the result of our classical calculation to the standard result obtained in the 
framework of fc^-factorization [4, 5]. 

It turns out that the formula we have just derived has a fairly simple con- 
nection with the result one would obtain in the framework of fc^-factorization. 

9 There are four bremsstrahlung diagrams of this type, and only one has been represented. 
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In order to see this connection, one needs the following relationship between the 
quantities p\ 2 (x±) an d the unintegrated gluon distribution 10 <pi l2 (k±): 

dtp(k ± ,x ± ) 2 V 2 (x±) 

<Pxl = 9 A ~k^~ ' (39) 

where dip(k±,x±)/d 2 x± is the number of gluons in a projectile, per unit of 
k\ and per unit area, and where d A = N 2 — 1 is the dimension of the adjoint 
representation. In fact, this identification of the unintegrated gluon distribution 
is more general than the model of eq. (34), for one has in general: 



dip(kj_,x±) 
d 2 x i 



J d 2 r ± e~ ik ^ (p a (x ± + r ± /2) Pa (x ± - r ± /2)) p , (40) 



where (• • -) p denotes the average over the hard color sources. From there, using 
eqs. (6) and (37), it is easy to write the pair production cross-section as: 



do\ 1 f d l ki± d 2 k 



!_ f d kl 

WJ W) 



*2± 



d<p 1 (ki±,x ± )di P2 (k 2± ,x ± -b) Tr (\ m ab + ( k i,k 2 ;q,p)\ 2 ^ 



I 



x / d z bd z x±- 



d 2 x±_ d 2 x±_ k 2 1± kl ± 



(41) 



where y q and y p are the rapidities of the quark and the antiquark. The inte- 
gration over b and x± is trivial to perform, and leads to the usual unintegrated 
gluon distributions: 

dai 1 f d 2 ki± d 2 k 



f d 2 k 1± d 2 k 2± . 

iww (li+21 " Pi_?i) 



dy p dy q d 2 p ± d 2 q ± (2wfd 2 A J (2tt) 2 (2tt) 

Tr ( | m ab + ( fc i^2;g,p)| 2 ) 

X^l(fel±)^(fe2±) * p—2 L ■ (42) 

This expression has a structure which is very similar to what one would have 
in the fcj_-factorization approach. Moreover, we have checked (see the appendix 

B) that the matrix element given by Tr ( |w~ b + (fci, k 2 ;q,p)\ ) is exactly iden- 
tical to the matrix element obtained in the framework of fc^-factorization by 
Collins and Ellis [4]. Note also that the x dependence here comes in through 
the unintegrated gluon distribution. This x dependence of unintegrated gluon 
distribution is related to that of the correlator of p's by eq. 40. The correlator 
in turn is determined by the weight function W in eq. 6 which, as discussed 
previously, satisfies a non-linear renormalization group equation in x. 

10 There are many definitions of the unintegrated gluon distribution in the literature. The 
one we are using here is: 



xG{x,Q 2 ) \ 



/ d(k 2 ± Mx,k ± ) . (38) 
Jo 
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From the previous formula, it is well known how to recover the standard 
results of collinear factorization. One must take the limit |fci_i_|, |fc2_i_| — ► in 
the quantity Tr ( |m~ b + (fci, /c 2 ; g,p)| ) /k\k\ (but not in the unintegrated gluon 
distributions). This limit is perfectly defined thanks to the Ward identities 
discussed in section 2.4. Then, the integration over the azimuthal angles of the 
vectors ki± and fc 2 ± gives the expression of the matrix element gg — > qq in the 
limit of collinear factorization, while the integration of the non integrated gluon 
distributions over kf ± and reconstructs the integrated gluon distributions. 
After this procedure, we are left with a factor S(p ± + q ± ) which naturally 
corresponds to the fact that in this limit the quark and antiquark must be 
produced back-to-back in the transverse plane. 



5 Beyond leading order 

Corrections due to terms of higher order in the hard color sources, breaking the 
k± factorization, are expected to become important in the soft regime, i.e. when 
the transverse mass of the quark or antiquark is of the order of the saturation 
momentum or smaller. In this section, we discuss the main issues that arise when 
one calculates corrections to the pair production amplitude that are of higher 
order in the hard color sources. In order to illustrate this discussion, let us first 
consider the example of corrections of order 0{p\pi) to the amplitude, which arc 
represented in figure 3. The contributions at this order can be grouped in three 





Figure 3: Contributions to the pair production amplitude at order 0{p\p2)- The 
diagrams involving bremsstrahlung have been omitted. The bold lines represent 
the hard color sources. 

sets, according to the number of classical field insertions along the quark line. 
The first set (the three diagrams in the upper left part of figure 3) involves only 
the Weizsacker- Williams fields A± and produced by the individual nuclei. 
The only difference with the previous order is that there are now two powers 
of and one power of A%. Then, there are two diagrams (represented in the 
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upper right part of figure 3) involving one power of the classical field A^ 2 and 
one power of the field A^. Finally, we have three terms with a single insertion 
of the classical field A^ 12 , i.e. the term of order 0{p\p?) in the solution of the 
classical Yang-Mills equation. 

Beyond leading order, there also arises the issue of the vacuum-vacuum di- 
agrams contained in the prefactor |(0i n |0 out )| 2 , since this factor cannot be re- 
placed by 1 anymore. Although this factor can be calculated order by order in 
an expansion in powers of the hard color sources, this makes the calculation of 
the single pair production probability much more difficult than the calculation of 
the average pair multiplicity. Alternatively, one could be less rigorous and sim- 
ply discard this prefactor in the calculation of P\ . This would of course violate 
unitarity (thereby resulting in a probability Pi larger than unity), which could 
then be restored in an approximate way by the replacement P\ — > P\ exp(— Pi). 
This way of restoring unitarity is equivalent to assuming that the distribution 
of pair multiplicities is a Poisson distribution. 

With this in mind, in order to perform the calculation of pair production 
beyond the leading order in the hard sources, one needs to go through the 
following steps: 

• solve the classical Yang-Mills equation with two sources up to the required 
order. This can in principle be done analytically by a generalization of the 
method presented in appendix A. But for a non-perturbative calculation 
to all orders, one has to resort to a numerical resolution of the classi- 
cal Yang-Mills equation with retarded boundary conditions. This step is 
under control after the work performed in [38, 39, 40, 41]. 

• compute the propagator of a quark in the previously obtained classical 
field. Only terms whose order is not higher than 0{p\p<i) need to be kept 
in the above example. For a systematic study of higher orders, one must 
rewrite the retarded amplitude u{q)T R v(p) in terms of retarded solutions 
of the Dirac equation, and then solve this equation numerically [42] . Here 
it is crucial to limit ourselves to the calculation of the pair multiplicity, 
which requires only the retarded quark propagator in the external field. 
Indeed, the numerical determination of the time-ordered propagator is 
a much harder problem, due to the great complexity of the boundary 
conditions obeyed by this type of propagator. 

• study the quantum evolution of the functional W[pi, p2]b]. Contrary to 
the leading order where only the correlators (pi. 2/01,2) were necessary, we 
need now higher order correlators. For instance, in the calculation at order 
0{p\p2) considered as an example, we would need the correlators (piPiPi) 
ad (pi pi pi pi). Indeed, strictly speaking quantum evolution does not pre- 
serve the Gaussian structure (unless in certain specific limits [46]). There 
are now ways to evaluate numerically these correlators [63] . Alternatively, 
at collision energies that are not too large, this step can be skipped since 
a Gaussian weight is a good approximation. 
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6 Conclusions 

We have discussed in this paper heavy quark pair production in heavy ion 
collisions, in the framework of the Color Glass Condensate model. We have 
computed explicitly the leading order contribution to the heavy quark pair pro- 
duction cross-section. This computation is valid in the kinematic region where 
fcj_ 3> Qs- We show explicitly that the classical result in this approximation 
is identical to the k± factorization results derived previously. We discuss how 
k± factorization may be violated for smaller k± as one enters the saturation 
regime. Numerical computations addressing this issue are in progress and will 
be reported separately [42]. 

We have not attempted here to relate our results to phenomenology. There 
have already been several phenomenological works applying k± factorization to 
various aspects of heavy quark production in Deeply Inelastic Scattering and 
in hadron-hadron scattering experiments. k± factorization ideas are also being 
studied in the framework of relativistic heavy ion collisions 11 . 
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A Classical color field to order 0(p\p2) 

We derive the classical field in the covariant gauge, and we follow closely the 
method of Kovchegov and Rischke. The gauge field is a solution of the classical 
Yang-Mills equation: 

[D„,F^] = r, (43) 
where J v is the classical current which must be covariantly conserved: 

[£>„, J"] = , (44) 

and the gauge condition is: 

9 M # = . (45) 

Making use of the gauge condition, the Yang-Mills equation can be rewritten 
as: 

UA V = J v + ig [A„ , + d^A v ] . (46) 

Let us first start with the terms and which are linear in the hard 
sources. Since the commutator [A^, F^ v + d^A v ] is at least quadratic in the 
sources, we can drop it at this order, and we have: 

□< 2 = ^2 , (47) 



^Private communication from D. Kharzeev and K. Tuchin. 
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where Jf 2 are the color currents associated to the individual nuclei. If we 
assume that the nucleus 1 is moving in the +z direction at the speed of light, 
while the nucleus 2 is moving in the — z direction, the currents are simply: 

J?, a = g6» + 6(x-) Plta (x ± ) , 

J^ a =g5^5{x + )p 2ta {x A _). (48) 

The index a is the color index carried by the classical sources and currents. This 
means in particular that at this order, the sources do not have any transverse 
components, i.e. that they describe recoilless objects. 

Eq. (47) then simply becomes a Poisson equation in the transverse plane, 
and the only non-zero components of A 1 ^ and A 2 can be written formally as: 

A t,a( x ) = -9$(x~)=L-piA x ±) > 

In fact, this solution will be needed in momentum space when we include it in 
the pair production amplitude: 

A+ a (k) = 2ng6(k-)-^Pi,a(k±) , 
k ± 

A 2 Jk) = 2ng6(k + )^ T p2, a (k±) ■ (50) 

In these formulas, k can be seen as the momentum flowing from the hard source 
to the quark line on which we insert the classical field. 

Let us now determine the contribution A^ 2 to the classical field. If we isolate 
the terms of Eq. (46) that are quadratic in the hard sources, we have: 

UA\ 2 = r i2 + ig [Api + A^Fr + + &>{A\ + A%)] , (51) 

where Jf 2 is the correction of order C(pip 2 ) to the color current. It can be de- 
termined by the current conservation which, including terms that are quadratic 
in the sources, reads: 

d v ji + d v r 2 + d v r 12 - 1 9 [a v1 + a v2 , j» + r 2 \ = o . (52) 

The current conservation at linear order in the sources tells us that d^J^ 2 = 0, 
so that we have simply: 

d v J v 12 = ig [A vl + A v2 , U{A vl + A v2 )] . (53) 

In addition to this equation, we know the following properties about Jf 2 : 

J^ 2 oc 8{x-)6{x+) , 
J12 oc S(x + )8(x~) , 

^2=0- (54) 
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The first two conditions mean that the + (resp. — ) component of the current 
has to move along with the nucleus going in the + (resp. — ) direction, and 
that effects of the second nucleus cannot start before the nuclei actually collide 
(hence the step functions). The third condition is simply a statement that the 
classical sources do not recoil. Eq. (53) can be made more explicit by writing 12 : 

d+ Jf 2 + d~J+ = ig [T 0) T b ] {A+ a DA- b - A~ b DA+ a } , (55) 

where T ajb are color matrices. Moreover, the current J± 2 , which follows the 
motion of nucleus 1, must depend locally on the source p\ and receive a non- 
local correction from the second nucleus. This means that in eq. (55) the term 
in A^ M DA 2b (which is non-local in p\ and local in p 2 ) must go into J^ 2 , and 
vice versa. At this point, we can write formally: 

J+=-ig[T a ,T b ]-^A- b nA+ a , 

J^=ig[T a ,T b ]^A+ a nA- b . (56) 

Note that the inverses l/d ± , which are a priori not uniquely defined, are made 
unambiguous by the step functions in eqs. (54). 

If we now realize that the only non-zero components of the first order 
strength tensor are: 

F[ + = -F+ l = d i Af , 

FT = -Fp = 5M 2 " , (57) 

the evolution equations for A± 2 read: 

UA+ 2 = J++ig[A 2 , d+A+] , 
DA^ = J u +ig[A+,d-A 2 ] , 

UA\ 2 = -ig [A+, 9M 2 ] - ig [A^^Af] ■ (58) 
Using the explicit form of the current Ji 2 , we can solve this as follows: 

A+ = - ig [T a ,T b ] 1 |(DA+ a ) ±A- b + (d+A+ a ) A",} , 

A u = ig[T a ,T b ] 1 1 (DA-,) ±A+ a + (d~A 2b "j A+ a | , 

A{ 2 = ig[T a ,T b }±{(d i AtjA- b -(d i A- b )Al a } . (59) 
Going to momentum space is trivial and leads to: 

AUk) = / ^ [ikt ~ §=) AiJk^A-^) , 

12 Wc denote d+ = d/dx' and d~ = d/dx+. 
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A\ 2 (k) = ig 1 ?^ J ^ {ik\ l k\) A+M^fo) , (60) 

where we use the shorthand k 2 = k — k\. One can note that the order 0(p\p 2 ) 
correction to the classical gauge field is proportional to a commutator of color 
matrices. This was to be expected as it is known that in an Abelian gauge 
theory there would not be such correction: the gauge field would simply be the 
sum of the gauge fields created by the individual nuclei. 



B Calculation of the traces 

In this appendix, we provide the result of the calculation of Tr Ijti^I 2 . In order 
to make this computation more compact, let us first introduce the following 4- 
vectors: 

a = (a + = 0, a~ =0, a± = q ± — feij_) , 
b= (b+ = 0,b-=0,b ± = k 1± -p ± ) , 



c + = p + + q^ 



k 2 
K l± 



k 



_,c = -p -q ,c± = k 2 ±~k 1± 

p +q p + + q + 



(61) 



and the following notations for the denominators: 

2q~p + + (q ± - k 1± ) 2 + m 2 = m 2 - i 



2q+p +(p ± -fci±) 
(p + qf = S , 



2 i 2 

1 m 



(62) 



where s, t and u are the standard Mandelstam variables for the gg — ► qq sub- 
process. We can rewrite the amplitude as follows: 



m ab + (k 1 ,k 2 ;q,p) = 
= ig 2 u{q) \ t a t b 



7 (m + 4h~* 



m 2 —t 



tbta 



7+(m + ^)7 i 
m 2 — u s 



>v(p) . 
(63) 



Depending on how terms are paired when squaring the amplitude, there are two 
kinds of color traces: 



tr c (t a t a t h t b ) = N C C 2 F 



tr c (t a t b t a tb) 



-lc. 



(64) 



where C F = (N 2 — 1)/{2N C ) is the Casimir in the fundamental representation 
of SU(N C ). In the large N c limit, the first trace scales like Af?, while the second 
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trace scales only like N c . We can therefore write the trace of the squared 
amplitude as: 



Tr (\m-+(k 1 ,k 2 ;q,p)\ 2 ) =g 4 C 
where we denote: 

7~(m + ^)7+ ^ 



N C C F T 3 - iri 



(65) 



T 3 = tr|(rf + m) 
+ (4 + m) 



m 2 -t 

7 + (m + J/)7~ i 
to 2 — u s 



+ - 

s 



7+(m + (jiYf 



m 2 -t 
7 _ (to + ^)7 + 



+ - 

s 



(66) 



and 



Ti = tr { (4 + to) 



7 (to + ^)7 ^ 
to 2 — £ s 



+ (^ + m) 



7 + (m + ^)7 ^ 



(jj/- m) 
— m) 



7 (m + ^)7 + 



7 + (m + (^)7 ^ 
to 2 - £ s 



(67) 



where the subscripts 1,3 refer to the order of the corresponding terms when 
N c — > oo. It is in fact a bit simpler to write 



Tr (|m-+(fci,fe;<jr,p)| 2 ) = ig 4 C F [N 2 C T 3 -T{] , 



(68) 



with T{ = Ti + T 3 , because it turns out that T[ has a more compact expression 
that T\. A direct calculation using form [64] leads after some rearrangement of 
the terms to the following expressions: 

, = 16 [ (m 2 - a 2 )p+q- | (to 2 - b 2 )p-q+ 



+ 



(to 2 - t) 2 (m 2 - u) 2 

(a-b-m 2 ){q ± -p ± -m 2 )-m 2 (qj_-p ± ) 2 + (p-a){q-b) + {p-b)(q-a)-2m 4 



(to 2 — <-)(to 2 — u) 



(69) 



and 



16(to 2 - a 2 )p+q- \&{m 2 - b 2 )p- q+ 16(q-c) 2 + 4c 2 s 
^ 3 — ; — — 1- 



(m 2 - t) 2 



(m 2 — u) 2 



+ 



+ 



s{m 2 - t) 



(a-c)(p + q -p q + + s/2) 



+ (2p + c -p ± -c±)(m 2 - q-a) - (2q c + - q ± -c^)(m 2 + p-a) 



s(m 2 — u) 



(b-c)(p+q- -p-q+-s/2) 
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+(2q + c - q ± -c±)(m 2 + p-b) - (2p c + - p ± -c±)(m 2 - q-b) . 

(70) 

It is straightforward to verify that T{ and T 3 vanish in the limits ki± — > and 
ki± -» p ± + q ± - 

One can also verify that our expressions, derived by solving perturbatively 
the classical Yang-Mills equation, are strictly equivalent to the formulas ob- 
tained in the framework of fcj_-factorization by Collins and Ellis [4] (also derived 
independently by Catani et al in [5]). In order to perform this comparison, one 
needs the following dictionnary (our notations are on the left, and theirs are on 
the right): 

q ^> p 3 , 

P Pi, (71) 

(the notation for the momenta of the incoming gluons is the same). The 
relationship between our amplitude squared and the formulas given in eqs. (5.7) 
and (5.9) of [4] is then summarized by: 

T[ = SY 2 , 

T 3 = 4(Y 1 +Y 2 ). (72) 

However, we do not agree with the amplitude squared given in the appendix of 
[7]. A private communication from Levin acknowledged the presence of a mis- 
take in the formulas quoted in [7] , but the corrected expressions were unavailable 
to us in order to compare them with our result. 
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